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Abstract 

The fluidic behavior of water at the micro/nano scale is studied by using of single DNA molecules 
as a model system. Stable curved DNA patterns with spans about one micron were generated by 
using of water flows, and observed by Atomic Force Microscopy. By rigorously comparing the 
numerical simulation results with these patterns, it is suggested that the form of the macroscopic 
hydrodynamic equation still works quantitatively well on the fluid flows at the nanoscale. The 
molecular effects, however, are still apparent that the effective viscosity of the adjacent water 
is considerably larger than its bulk value. Our observation is also helpful to understand of the 
dynamics of biomolecules in solutions from nanoscale to microscale. 
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Novel micro-fluidics has provided new possibilities for the development of fabricating 
high-performance devices with micrometer or even submicrometer dimensions , as 

well as bio-processing for bio-medical analysis 0. The ability to control accurately the 
flow requires the understanding of fluid flows at the micro/nano scale near macro molecules 
or above solid surfaces. Recent experiments have shown that macroscopic hydrodynamics 
can reliably describe the flow of fluids through channels with cross-sectional dimensions 
that range from tens to hundreds of micrometers 51. The nanometer scale, however, is in 
the transition regil „e between continue and tnolecnle donated conditions QQ, aoont 
which we know very little 6\. Molecular dynamics has been used to study the nanofluid 
[7] and consistence between the numerical results from molecular dynamics and Navier- 



Stokes equations has been obtained recently 



However, both high and normal 'effective' 



viscosities, its departure from the bulk value usually characterizing the molecular effects, 
and boundary slip have been observed in the confined water when the sample thickness is 
comparable to or below the diameter of DNA molecules [Oj, indicating that the dynamics 
of nanofluid should be much complex. Moreover, it has been found experimentally that the 
water flow th roug h orifices with diameter about 8 /im is considerably different from that of 
ordinary size [lfj. At a hydrophobic surface, in an unexpectedly extended (~ 4 nm) water 
aver, the density of water has a noticeable reduction (85-90% of the density of bulk water) 



11 1 . Consequently, the fluidity of water several to tens of nanometers above non-liquid 



materials should be different from that of the macroscopic systems despite the fact that the 
water molecules are much smaller. 

The difficulty in the experimental studies is partially due to the limited experimental tech- 
nologies to obtain data that can be compared by theoretical predictions quantitatively. Al- 
though molecules with fluorescent labels have been widely used [3] to make the biomolecules 
visible, the resolution is limited by the minimal photon wavelength, which is about half a mi- 
cron. Atomic Force Microscopy (AFM) and Electronic Microscopy (EM) have the nanometer 
resolution. However, they require the objects to be fixed on solid substrates. It is obviously 
questionable that the quantitative details on the dynamics can be exploited from motionless 
objects unless they are in steady-states, which are difficult to be confirmed. 

The study of biological systems will continue to inspire the development of new physics 



13 1 and the dynamics of DNA in fluid flow has been extensively studied by various methods 
12.[jjJ. Thanks to the elastic behavior of DNA fibers that they reach maximal lengthes with 
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a tension larger than ~65 pN j3| and are not broken until 400 pN Q| , those overstretched 
DNA fibers can reach steady states in the fluid flows and be fixed on substrates. In this 
Letter a simple experiment was then set up to obtain steady curved DNA patterns for 
span lengths about or less than 1 /im, generated by water flow. By rigorously comparing 
the observed DNA patterns with theoretical predictions, it is suggested that macroscopic 
hydrodynamic equations still works quantitatively well on the fluid flow at the nanoscale at 
least at steady states. To our best knowledge, this is the first report on the applicability of 
the Navier-Stokes equations on the fluidic behavior of water near single DNA molecules at 
the nanoscale by testing the numerical predictions with the experimental observation. 

The ADNA, purchased from Sigma Co. (USA), was diluted before use to a concentration 
of 3 ng//il with TE buffer (40 mMol/L Tris-HCl, 1 mMol/L EDTA, pH 8.0). The experiment 
jegan by depositing a drop of 10 /d ADNA solution onto a mica sheet modified by APTES 



171 ] that was placed on a rotary device. The distance between the drops of ADNA solution 
and the axis of the rotary device was always about 10 mm. The spin speed of rotation is 
kept to be 1150 rpm. After rotation, the drop would be thrown off and some DNA molecules 
were aligned on the surface. To demonstrate that DNA molecules were aligned by water flow 
rather than the liquid-gas interface, some DNA fibers were observed by Tapping- Mode AFM 
(Nanoscope Ilia, Veeco Instruments, Santa Barbara, CA) using silicon tips (Silicon-MDT 
Ltd., Russia) jj. It is found that the relative extensions of those aligned DNA molecules 
are within 100% ± 5% while those are about 133% to 150% in the molecular combing by the 
receding meniscus Then another drop of water was placed on the same area of the 

surface on which there are aligned DNA molecules while the orientation of the mica sheet 
on the rotary device was changed about 90°. A subsequent rotation created a second flow 
on the mica surface, approximately perpendicular to the first one. The velocity of the drop 
was estimated to be about 1 cm/s, with an upper limit of 3 cm/s. The as-generated DNA 
patterns were revealed by Tapping-mode AFM jjj. 

We have systematically observed over 100 curves with span length L, 50 nm < L < 2 
/im after the second water flow, as typically shown in Fig. 1. Some of the DNA segments 
still have their original aligned patterns, while other segments bend in the direction of the 
second water flow. Among them, 7 symmetric curves are constituted of those overstretched 
DNA fibers, which reach their maximal lengthes, with span length ranging from ~580 nm 
to ~1260 nm. The fact that the ratio of the contour and the span is about 1.75 for those 
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FIG. 1: DNA patterns obtained by AFM imaging. Some of the DNA segments still keep their 
original aligned patterns, while some segments bend to the moving direction of liquid flow and 
form wave-like patterns. The curves denoted by red dots are the theoretical predictions 5ft from our 
numerical calculation which agree extremely well with the curves from the experiment. The light 
magenta arrows are only served as a guide to eyes rather than the exact directions of the second 
water flows. The blue dots in the inset are the data extracted from the DNA pattern and the red 
line is the theoretical prediction 3?. 

overstretched patterns further verifies that DNA molecules are aligned and curved by water 
flows rather than the interface. The relative extension is calculated from the contour and 
the span length L of a curve. In the following we will demonstrate that these symmetric and 
steady DNA curves are a result of the balance of the macroscopic hydrodynamic forces with 
an adjusted viscosity, as illustrated by the 5 typical examples shown in Fig. 1, i.e., curve 2 
in (a) and curves 1, 2, 3, 5 in (b). The other patterns, including the asymmetric curve 1 in 
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Fig. 1(a) and the double-peaked curve 4 in Fig. 1(b), have been initially analyzed in 
and will be discussed in details elsewhere. 

The numerical analysis is based on the bead-spring model [^J. To test whether the 
macroscopic hydrodynamic equations still works quantitatively well, the forces acting on 



the beads are computed by simulating water flow according to the Navier-Stokes equations. 



Explicitly, the lattice Boltzmann simulation 



2l| is applied to simulate the water flow around 



the spherical beads. A non-slip boundary condition [22| is applied at the bead-fluid bound 



aries. The momentum exchange 



23, 



24j at the solid-fluid interface is used to determine the 



force acting on the DNA beads. In the lattice Boltzmann simulations, the non-dimensional 
radius of each bead is 5.2 lattice units and the relaxation time r = 0.75. The density and 
kinetic viscosity of the water are 998.3kg /m 3 and 1.007xl0 _6 m 2 /s j^, respectively. 




FIG. 2: The system for our numerical simulation together with the coordinates. The red circles 
are the beads placed on the DNA pattern scaled down from the curve marked by 2 in Fig. 1(a). 
The black arrows represent the directions and magnitudes of the hydrodynamic forces on the beads 
obtained by numerical simulations. There are three forces acting on the DNA segment between C 
(x = 0) and x = 5x: the tension Tq at C along the x-axis, the tension Ti at the other end of 
the segment and the hydrodynamic force fh{5x) on this segment, which can be computed from the 
hydrodynamic force distribution. 



Let us take the curve marked by 2 in Fig. 1(a) as an example. The settings of our 
numerical simulations are shown in Fig. 2. Some beads, say sixty-nine, were placed on the 
DNA pattern scaled down from the curve. The radius of each bead is 1 nm, approximately 
equal to that of the DNA strands, while the center-to-center distance between the nearest- 
neighbor beads is always 2.308 nm. Periodic boundary conditions are applied for the x and 
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z directions. The mesh sizes are: L x is a little larger than the span length L of the DNA 
pattern, the inlet and outlet in y direction is always 30.5 nm from the closest bead, see, e.g., 
L x = 95.5 nm and L y = 120.0 nm for L=90.9 nm; L 2 =41.1 nm. It is found that L z is large 
enough so that the effect of the periodic boundary condition in the z direction is negligible. 
The fluid velocity at the inlet and outlet is assumed to be 1 cm/s and along the -y direction 
at all time. 

The force distribution acting on beads is shown in Fig. 2. Since the actual DNA fiber is 
continuous rather than discrete, we therefore used the best continuum curve fitting from the 
calculated discrete force to approximate the realistic hydrodynamic force distribution along 
the DNA fiber. 

It is evident that if a DNA fiber is in a steady state of the hydrodynamic force, the total 
force, i.e., the sum of the hydrodynamic force and the tensions, acting on any segment of 
the DNA fiber should vanish. Unfortunately, the tension is unknown and it is difficult to 
measure the velocity of water flow acting on the DNA molecules to an acceptable accuracy 
so that we cannot examine the equilibrium condition explicitly. However, if we assume that 
the coordinates of the two ends of the DNA at locations C and O (see Fig. 2) are given, 
then the hydrodynamic force distribution obtained above will uniquely determine a curve 3ft. 
Further, 3ft can be used to compare with the experimental result. If these two curves coincide 
with each other, the DNA pattern obtained by the experiment is indeed in the steady state 
of the hydrodynamic force. Otherwise the discrepancy reflects the special behavior at the 
nanoscale. 

To construct the curve 3ft, we start our steady force balance calculation from the critical 
point C as shown in Figure 2. Assuming that there is a small displacement Sx in the re- 
direction, then there are three forces acting on the DNA segment between C (x = 0) and 
x = Sx: the tension at C along the x-direction, the tension T x at the other end of the 
segment and the hydrodynamic force ih{Sx) on this segment. 8x = 0.5 nm was used in our 
calculations. For the equilibrium state, the force balance for the DNA segment gives: 

T c + Tx + i h (Sx) = 0. (1) 

If we assume that the magnitude of |Tc| is Tc, both the magnitude and direction of the 
tension Ti can be obtained through Eq. [T] Consequently we can obtain the y coordinate 
of the DNA segment at x — Sx. Next, we consider the position of the curve 3ft at x = 25x. 
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The tangential direction at this point can be obtained similarly to find the increment of 
y value from x = 5x to x = 25x. This extrapolation process can be repeated until x a is 
reached. Then an iterative process is used to self-consistently determine Tc to meet the 
value y required at O. In this way, we obtain the curve 3ft. It should be noted that the force 
to bend the DNA molecules to the curves shown in Fig. 1, except for the parts near the 
anchored points, is about 0.1 pN, negligibly smaller than the tension in the DNA curves. 

3ft is displayed as red circles in Fig. 1(a). To our great surprise it is found that 3ft from our 
theoretical prediction agrees extremely well with the curve from the experiment. Numerical 
simulations for 35, 49 and 89 beads, corresponding to span lengths L = 45.1, 63.6, and 116.6 
nm, respectively, have also been carried out. It is found that the normalized patterns from 
numerical predictions with different length scales can be normalized to a high accuracy. In 
the same way, we obtain the theoretical predictions 3ft for the curves 1, 2, 3, 5 in Fig. 1(b). 
In order to obtain a quantitative comparison between the numerical predictions and the 
experimental results, the global relative error is computed from those 5 curves, which is 5 
± 1 %. 

The tension Tc at the critical point C from our calculations for different L is shown in 
Fig. El which fits very well to a linear function Tc = 0.0156L when the inlet fluid velocity 
equals 1 cm/s. Following this linear scaling, we found that Tc = 9.0 and 19.7 pN for L = 580 
and 1260 nm, respectively. 

The tensions are larger than 65 pN in the overstretched DNA fibers with relative 
extension about 1.75. However, Tc calculated for L = 580 nm is only 9.0 pN. Even for the 
upper limit of the velocity, 3 cm/s, of the water drop, Tc = 27.0 pN is still smaller than 65 
pN. We also note that the viscosity of water is the same as its bulk value in the calculation. 
If we assume that the local effective viscosity near the DNA fiber is more than 3 times of its 
bulk value, the calculated tension in all the 7 steady and symmetric curved DNA patterns 
will be larger than 65 pN. Enhanced by a factor of 6 for the effective viscosity of water near 
proteins has been also observed recently j^. In the inset of Fig. El we have included the 
tension distribution along the DNA fiber for L of 1260 nm. 

The numerical analysis in the manuscript has neglected the fact that DNA fibers are close 
to the mica surface. It has shown that polymer molecules can move on the surface with a 
low barrier although the absorption on molecules by the surface is strong [27]. Therefore, 
we believe the interaction between the extended DNA strands and the substrate will not 
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FIG. 3: The tension Tc at the point C of the DNA curve shown in Figure 2 for different span 
lengths, and their linear fit function. The inset is the tension distribution along the DNA fiber for 
the span length of 1260 nm, computed from a bulk viscosity. The tension at C has its minimum 
value and the tension at O has its maximum value To ~ 1.9Tc. 

change the conclusion we obtain above. In this case the flow around the DNA fibers may be 
a shear flow. We have performed new numerical simulations. The distance from the x — y 
plane, where the DNA curve located, to a non-slip boundary is (. The velocity of the flow at 
the inlet increases linearly from the boundary whereas the velocity in the line with the x — y 
plane is 1 cm/s. The other settings are the same as before. It is found that the normalized 
patterns 3? obtained in this way for ( = 3.84 nm and 0.58 nm are consistent with those 
obtained above perfectly, while the tensions increase by a factor of 1.2 and 2.6 separately. 
The fluid velocity near wall is usually much smaller than that of the drop of water, therefore 
higher effective viscosity is also demanded. 

The accuracy of the fitting of the theoretical predictions to the experimental patterns is 
at the level of AFM resolution, about 10-20 nm. This observation suggests that the force 
exerted by the adjacent water on the DNA fiber can be determined quantitatively by the 
macroscopic hydrodynamic equations at least in the steady state even at a resolution of 
about 20 nm. The viscosity of the adjacent water, however, is much larger than its bulk 



S 



value. It should be noted that the excellent agreement between our numerical predictions 
and the experimental results further demonstrate the robustness of the numerical model 
although the detail of the experiment should be much complex. 

We thank Dr. Gary D. Doolen, Profs. D. Bensimon, V. Croquette and Z. Ye for helpful 
discussion and suggestions, Miss Chunmei Wang for extracting DNA data from AFM images. 
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